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ABSTRACT 

The deuterium abundance evaluation in the direction of Capella has for a long time been used as a 
reference for the local interstellar medium (ISM) within our Galaxy. We show here that broad and weak 
Hi components could be present on the Capella line of sight, leading to a large new additional systematic 
uncertainty on the JV(Hi) evaluation. 

The D/H ratio toward Capella is found to be equal to 1.67 (±0.3) x 10~ 5 with almost identical \ 2 
for all the fits (this range includes only the systematic error ; the 2 a statistical one is almost negligible 
in comparison). It is concluded that D/H evaluations over Hi column densities below 10 19 cm~ 2 (even 
perhaps below 10 20 cm~ 2 if demonstrated by additional observations) may present larger uncertainties 
than previously anticipated. It is mentionned that the D/O ratio might be a better tracer for Di 
variations in the ISM as recently measured by the Far Ultraviolet Spectroscopic Explorer (FUSE). 

Subject headings: cosmology: observations — ISM: abundances — Galaxy: abundances — 
Ultraviolet: ISM — stars: individual (Capella) 



1. INTRODUCTION 

Deuterium is understood to be only produced in signif- 
icant amount during primordial Big Bang nucleosynthesis 
(BBN) and thoroughly destroyed in stellar interiors. Deu- 
terium is thus a key element in cosmology and in galactic 
chemical evolution (see e.g. Audouze & Tinsley 1976). In- 
deed, its primordial abundance is the best tracer of the 
baryonic density parameter of the Universe f2s, and the 
decrease of its abundance during the galactic evolution 
should trace the amount of star formation (among other 
astrophysical interests). 

In the Galactic ISM, D/H measurements made toward 
hot stars have suggested variations : IMAPS observations 
toward S Ori led to a low value (Jenkins et al. 1999), con- 
firming the previous analysis by Laurent et al. (1979) from 
Copernicus observations, while toward 7 2 Vel they led to a 
high value (Sonneborn et al. 2000). This seems to indicate 
that in the ISM, within few hundred parsecs, D/H may 
vary by more than a factor ~ 3. 

In the nearby ISM, the case of G191-B2B was studied 
in detail (see the most recent analysis by Lemoine et al. 
2002) and the evaluation toward Capella (Linsky et al. 
1995) taken as a reference. Their comparison provided, 
for a while, a possible case for D/H variations within the 
local ISM. 

Concerning G191-B2B, Lemoine et al. (2002) have 
shown that the total A to t(Hi) column density evaluation 
was greatly perturbed by the possible addition of two 
broad and weak H I components. Such components, able 
to mimic the shape of the Lyman a damping wings, can in- 
duce an important decrease of the evaluated A to t(Hi). To 
illustrate this point, the error bar estimation on AT tot (Hl) 
from all previously published studies considered as the ex- 
tremes of a 2a limit was of the order of dcx 0.07, while 
including the Lemoine et al. (2002) analysis enlarged the 
error bar to about dex 0.37. This huge change has, of 



course, a considerable impact on any D/H evaluation. 

This raises two crucial questions. First, is that situation 
typical of G191-B2B alone and possibly due to an unex- 
pected shape of the core of the stellar Lyman a profile 
improperly described by the theoretical models? Second, 
if weak H I features are present in the ISM, to what extent 
are evaluations toward other targets affected? 

2. SUMMARY OF THE G191-B2B CASE 

From the combination of STIS echelle observations 
(spectrograph on board the Hubble Space Telescope, HST) 
and FUSE ones (the Far Ultraviolet Spectroscopic Ex- 
plorer, Moos et al., 2000), Lemoine et al. (2002) have 
found through iterative fitting process (with the Owens . f 
fitting program developed by Martin Lemoine and the 
French FUSE team) that three interstellar absorption com- 
ponents are present along the line of sight and that two 
additional broad and weak H I components could be added, 
detected only over the Lyman a line (negligible over the 
Lyman (3 line) but important enough to strongly perturb 
the total H I column density evaluation. 

Within the local ISM, it has been shown that such addi- 
tional Hi absorptions are often present ; they have been in- 
terpreted either as cloud interfaces with the hot gas within 
the local ISM (Bertin et al 1995) or as "hydrogen walls" , 
signature of the shock interaction between the solar wind 
(or stellar wind) and the surrounding ISM (Linsky, 1998). 
This latter heliospheric absorption has been modeled by 
Wood et al. (2000) and a prediction derived in the direc- 
tion of G191-B2B (see Figure 9 of Lemoine et al. 2002). 
Most of the predicted absorption is expected in the satu- 
rated core of the observed interstellar line but some weak 
absorption (~ 5% of the continuum) might extend over 
several tenths of angstroms on the red side of the line, due 
to the neutral hydrogen atoms seen behind the shock in 
the downwind direction where G191-B2B is located. It 
was found that the combination of two broad and weak 
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Hi components can easily reproduce the model prediction. 
If real, besides the three interstellar absorptions, a fourth 
component representing the bulk of the predicted absorp- 
tion and a fifth one for the broad and shallow extended 
red wing are needed. This is exactly what Lemoine et al. 
(2002) have found. 

In the course of determining the minimum number of 
components (each defined by its Hi column density N, its 
velocity v, its temperature T and turbulence broadening 
£) needed to fit the data, Lemoine et al. (2002) completed 
the F— test which uses the Fisher-Snedecor law describing 
the probability distribution of \ 2 ratio. What is tested is 
the probability that the decrease of the \ 2 with additional 
components is not simply due to the increase of free param- 
eters. The result gives a probability < 10~ 4 and ~ 5 10 -3 
that a fourth and a fifth Hi component are respectively not 
required by the data. These low probabilities of non oc- 
curence strongly suggest that Lemoine et al. (2002) have 
indeed detected the heliosphcric absorption downwind in 
the direction of G191-B2B. 

Note however that this heliospheric complex absorption 
profile is simulated by two components whose physical 
meaning in terms of hydrogen content and/or tempera- 
ture is not clear. Furthermore, the photospheric Lyman a 
stellar core is difficult to evaluate (see discussion in e.g. 
Lemoine et al. 2002) and is slightly red-shifted relative to 
the ISM absorptions ; this result may very well be simply 
related to the use of a white dwarf as background target 
star. The detailed analysis of the Capella line of sight 
could directly test the heliospheric hypothesis. 

3. THE CASE OF CAPELLA 

If the two additional components present along the 
G191-B2B line of sight are as a matter of fact due to 
an heliospheric phenomenon, it is an extremely local sig- 
nature (within few hundreds of astronomical units to be 
compared to the few tens of parsecs lines of sight lengths) 
which should be also present along the Capella sight-line, 
both stars being separated by only 7° on the sky, and sim- 
ilar in shape to the structure predicted and observed in 
the direction of G191-B2B. If that description is correct, 
we are expecting an extra absorption reasonably repre- 
sented by two additional components, a main one mostly 
lost within the ISM absorption core and a weak one ex- 
tending over several tenths of angstroms on the red side 
of the line, again due to the neutral hydrogen atoms seen 
behind the shock in the downwind direction where both 
G191-B2B and Capella are located. 

Recently, Young et al. (2002) analysed new obervations 
obtained at Lyman /3, Lyman 7 and the whole Lyman 
series with FUSE. The precise Lyman a stellar profile 
compatible with all Lyman lines and with the data sets 
obtained at different phases of the Capella binary system 
(see also Linsky et al. 1995) was reevaluated (Wood, 2001) 
and is used here as a reference profile, Sr. 

We thus revisited the fits completed over the Lyman a 
line as observed toward Capella with the best available 
data set, i.e. the one obtained with the GHRS (the God- 
dard High Resolution Spectrograph on board HST). The 
study by Linsky et al. (1995), essentially confirmed by 
Vidal-Madjar et al. (1998), shows that only one interstel- 
lar component (the Local Interstellar Cloud, LIC, also seen 



toward G191-B2B) is needed on that line of sight. This 
very simple structure strengthens the Capella case as the 
simplest one where D/H can be very well evaluated. How- 
ever, Vidal-Madjar et al. (1998) have already noted that 
an additional weak and broad Hi component was required 
to better reproduce the profile ; this was a first indica- 
tion of the presence of an heliosphcric absorption toward 
Capella. In fact, we were able to show that, as in the 
case of G191-B2B, the addition of one or two weak and 
broad Hi components (together with the very weak geo- 
coronal component present at a known velocity but not 
shown on Figure 1 for clarity) improves the % 2 . More pre- 
cisely we fitted the GHRS data assuming that the stellar 
continuum was Sr. Adding successively to the fit one then 
two free Hi components (the added Hi components have 
only three free parameters, velocity v, column density N 
and width T, since the thermal T or turbulent broaden- 
ing £ act in an undifferentiated manner when only one 
species is observed) we obtained the following \ 2 /degree 
of freedom (d.o.f.) values: for only the LIC component and 
the geocorona, 844.89/716 ; for one additional component, 
831.17/713 ; for two additional components, 822.68/710. 
The F— test probabilities that these two additional compo- 
nents are not required by the data are respectively 8.5 10~ 3 
and 6.3 10 -2 . The first one is clearly needed here (its cor- 
related parameter ranges according to different possible so- 
lutions similar in terms of x 2 are : 12 < v (km s _1 ) < 22 ; 

1.4 10 14 < N (cm- 2 ) < 4.4 10 15 ; 10000 < T (K) < 32000) 
but unlike in the case of G191-B2B, the second one corre- 
sponding to the weaker and broader one (parameter ranges 
are : 24 < v (km s" 1 ) < 40 ; 3.6 10 12 < N (cm~ 2 ) < 

8.5 10 13 ; 120000 < T (K) < 300000) is less strongly 
needed. These ranges are certainly compatible with the 
corresponding estimated values in the direction of G191- 
B2B (see Figure 10 of Lemoine et al. 2002). 

To search for the possible impact of the choice of the 
continuum on the evaluation of iV(Hi), we fixed this value 
and looked for the best fitted solutions while the stel- 
lar continuum we used Sr was allowed for some varia- 
tions by multiplying it by a low order polynomial (8 th or- 
der) which coefficients were free to vary along with all 
components parameters. Results are shown in Figure 1. 
Slight changes of the continuum shape by no more than 
±10% lead to nearly identical \ 2 values, with AT(Hi) vary- 
ing from 2.0 x 10 18 cm~ 2 to 1.4 x 10 18 cm~ 2 which corre- 
sponds to a change in D/H from 1.37 x 10~ 5 to 1.96 x 10~ 5 . 
This is clearly a larger range (~ ±0.3 x 10~ 5 ) than the one 
previously claimed (< ±0.2 x 10~ 5 , Linsky et al. 1995). 

The situation could be even worse since we do not 
know how far the Capella continuum could be away from 
Sr. The question is thus to evaluate if the Capella Ly- 
man a stellar continuum shape is estimated to better than 
±10% or not ? It is true that having a binary system can 
help constraining the continuum shape as Linsky et al. 
(1995) did, but their whole approach requires that the Ly- 
man a stellar profiles of both Gl III and G8 III stars are 
invariant with phase and time. In fact, from the study of 
120 IUE echelle spectra, Ayres et al. (1993) have shown 
on one hand that the line fluxes were surprisingly stable, 
but on the other hand that whichever way they process the 
data, obvious variations were seen. These seem to be re- 
lated to variations of the blue peak of the Gl III dominant 
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stellar Lyman a line. They found that in the 1981-1986 in- 
terval, the line shape at phase 0.25 of the system was quite 
stable and similar to the one recorded with the GHRS in 
1991 (at a ±10% level). Earlier spectra taken in 1980 or 
later ones observed after 1986 look quite different. This 
very careful study shows that with the IUE sensitivity level 
of ±5%, variations are clearly detected. Since Linsky et al. 
(1995) used GHRS observations at two different phases of 
the system (0.26 and 0.80) taken respectively in April 1991 
and in September 1993, i.e. two and a half year apart, it 
is difficult to ascertain that the Lyman a profile evaluated 
for each stellar component is well controlled. Because of 
the very careful analysis made by Linsky et al. (1995), 
it may be possible that the stellar Lyman a profiles are 
relatively well evaluated but certainly not at a level better 
than ±10% as previously mentioned. 

Thus, an hcliosphcric absorption is also detected on the 
Capella sight-line ; furthermore, even in such a simple 
ISM configuration (a unique component), it appears im- 
possible to tightly constrain the total Hi column density 
in that direction. 

4. DISCUSSION 

We have shown that, for two lines of sight, AT (Hi) can- 
not be evaluated with a high accuracy. Column densities 
on both sight-lines are very similar, of the order of few 
times 10 18 cm . For lower column densities, the situ- 
ation should be worse since then the possible absorption 
signature of the weak components is becoming relatively 
more and more important and the Lyman a line is getting 
closer to the flat part of the curve of growth where column 
densities are indeed difficult to evaluate. 

Note however that the HST/EUVE comparison com- 
pleted by Linsky et al. (2000) shows that often AT (Hi) 
values derived from the GHRS and EUVE data (not sen- 
sitive to weak Hi) are in good agreement, implying that 
hcliosphcric absorption (or other hot components) don't 
necessarily ruin Lyman-alpha analyses in a dramatic way. 
But clearly counter examples leave that question open due 
to possible systematics related to the evaluation of total 
Hi below the Lyman limit in the EUVE domain. 

On the contrary, one could guess that for larger column 
densities the situation should improve since the Lyman a 
damping wings are becoming broader and the signature of 
the weak features may disapear in the line core. Just above 
10 19 cm -2 , the reliability of the D/H values is greatly en- 
hanced if the studied gas is demonstrably warm (6000 K ; 
for a thorough discussion see York, 2001) ; as one goes 
above 10 20 cm -2 , credibility increases unless either cold 
gas components are hidden in the warm Di but still affect 
the Hi damping wings or weak Hi features at high velocity 
are present. 

The unknown referee further stressed this point through 
an impressive report. He mentionned that at N(Hi)~2. x 
10 19 cm -2 , the half-intensity point of pure damping Ly- 
man a and Lyman (3 profiles, located respectively at ve- 
locity shifts of 274 km s _1 and 55 km s _1 , should be the 
place where a putative high-velocity feature could have a 
strong perturbing influence on the damping profile. How- 
ever only very few high velocity ISM components were 
detected above 120 km s _1 . This could lead to the in- 
verse impression that for larger column densities, the es- 



timation through the Lyman (3 line should be more ques- 
tionable than the one made at Lyman a. As a matter of 
fact, in both the 7 Cas and ( Pup lines of sight, the Hi 
column density was discrepant when derived through the 
Lyman a or the Lyman (3 line (see below). However in 
both cases the Lyman (3 estimations of AT(Hi) are smaller, 
in contradiction with the formerly suggested cause since 
the most perturbed evaluation by additional absorptions 
should lead instead to larger column densities. 

High velocity ISM components essentially observed be- 
low 120 km s _1 were only searched for through other lines 
and species than Hi at Lyman a the strongest transition 
of the most abundant element. For instance, Cowie et 
al. (1979) reported Lyman S and Lyman e Hi ISM absorp- 
tions up to about 105 km/s for t Ori. From their study, the 
referee evaluated that a shock at 274 km s _1 should pro- 
duce either a very broad (~ 175 km s _1 ) and undetectable 
(maximum depth of 5.9 x 10 -3 ) post-shock absorption sig- 
nature or should originate in a region far downstream from 
the front where the gas has cooled and compressed enough 
to allow recombination of the H atoms, i.e. a shock from a 
supernova explosion entering the radiative phase. In this 
second case however, he estimated from Cowie and York 
(1978) and Spitzer (1978) that such signatures should oc- 
cur only very near a known supernova event, i.e. within 
about 30 to 60 pc for standard ISM and SN values. Thus, 
that looks unlikely too. 

On the other hand, high- velocity gas could be generated 
by the target stars. While Gry, Lamers and Vidal-Madjar 
(1984) seem to detect most of the activity at velocities 
below 100 km/s, they nevertheless identified a transient 
component at -150 km/s toward j 2 Vel through Lyman 5 
and another one at -220 km/s toward t Ori, through Ly- 
man e. Note also that this survey was completed over 
a limited spectral domain scanned with the Copernicus 
instrument, and not at Lyman a, i.e. with a relatively 
limited sensitivity. 

One might argue that there are some risks that stellar 
ejecta could influence the Lyman a measurements ; but 
the observers have a good defense : multiple observations 
at very different epochs. This strategy was invoked by 
Jenkins et al. (1999) and Sonneborn et al. (2000) in their 
studies of D/H toward 5 Ori, £ Pup and -f 2 Vel. Their 
findings are thus pretty convincing in this regard. 

All the above stated arguments should mitigate our con- 
cern that small amounts of Hi at high velocities are a likely 
source of confusion for the flanks of the damping profiles 
for A 7 " (Hi) of the order of or greater than 2. x 10 19 cm~ 2 . 

One however should recall the two lines of sight for 
which the Hi column density was discrepant when derived 
through the Lyman a or the Lyman f3 line : 

• 7 Cas 

in Bohlin, Savage and Drake (1978), Lyman a only 
(Copernicus) AT(Hi) = 1.45 (±0.29) x 10 20 cm" 2 ' 

in Ferlet et al. (1980), core of Lyman f3 only (Coper- 
nicus) AT(Hi) = 1.10 (±0.10) x 10 20 cm- 2 ' 

in Diplas and Savage (1994), from Lyman a only 
(IUE) AT(Hi) = 1.48 (±0.31) x 10 20 cm" 2 

• C Pup 
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in Bohlin (1975) from Lyman a only (Copernicus) 
iV(Hl) = 0.97 (±0.05) x 10 20 cm" 2 
in Vidal-Madjar et al. (1977), from Lyman f} only 
(Copernicus) A(Hl) = 0.80 (t° \) x 10 20 cm" 2 
in Diplas and Savage (1994), from Lyman a only 
(IUE) JV(Hl) = 0.89 (±0.17) x 10 20 cm" 2 

in Sonneborn et al. (2000), from Lyman a only 
(IUE) JV(Hl) = 0.92 (±0.01) x 10 20 cm" 2 

Since Lyman /3 is less sensitive to weak Hi features, it is 
interesting to note that in both cases the Lyman [3 eval- 
uation is slightly below the Lyman a one, a possible in- 
dication of a similar effect for column densities of the or- 
der of 10 20 cm -2 . Note however that these two cases are 
marginally convincing since the different evaluations are 
still compatible within the error bars. 

Therefore, if this effect is indeed real for relatively large 
column densities, it means that some of the D/H values 
may be underestimated and thus the higher D/H ratios 
may be favoured. This further shows the importance of 
the 7 2 Vel estimation (Sonneborn et al. 2000). 

5. CONCLUSION 

According to our conclusion that there is lower precision 
for A(Hi) measurements within the local ISM, i.e. at low 
column densities, which induces large error bars on the 
corresponding D/H estimations: 

• the A (Hi) evaluation in the direction of Capella is 
relatively less accurate than previously claimed and 
is of the order of log TV (Hi) = 18.22 (±0.08) leading 
to D/H = 1.67 (±0.3) x 10~ 5 ; 

• an average D/H ratio may exist in the local ISM but 
should be larger than previously evaluated since lo- 
cally iV(Hi) could be overestimated by as much as 



about 20% ; this does affect arguments about local 
variability ; 

• D/O should be a better tracer of D variations as 
originally suggested by Timmes et al. (1997), and 
directly verified in the LISM (Moos et al. 2002 ; 
Hebrard et al. 2002a, 2002b) and further confirmed 
by the stability of O/H over longer path length 
(Meyer et al. 1998 ; Andre et al. 2002), provided 
all sources of errors related to the oxygen measure- 
ments are well understood. 

Finally, we note also that a similar systematic effect has 
been pointed out by Pettini and Bowen (2001) for the eval- 
uation of D/H in quasars absorption line systems. Again 
as in our study, the systems presenting the highest values 
of A(Hi) are derived from the damping wings of the Ly- 
man a line, which also includes all Hi in close proximity 
of the hydrogen at the line center of deuterium (Buries, 
2001), while those for which A (Hi) is evaluated from the 
discontinuity at the Lyman limit present smaller column 
densities and larger D/H evaluations. 

Clearly, many more lines of sights should be analysed to 
resolve this issue. 
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Fig. 1. — GHRS observations of the Capclla Lyman a line (thin histogram). Since Capclla is close to G191-B2B on the sky, two broad and 
weak Hi absorption components (thick dashed lines) were freely added in the fitting process. The outcome is very similar to the one found 
for the G191-B2B line of sight. In each panncl, from top to bottom, the Hi column density of the unique ISM component is fixed at different 
values (noted on the upper right corners). The fits are shown as a thin almost invisible line running through the data points (according to 
the very high quality of these fits) in all three pannels. Since these fits are very similar in terms of \ 2 it is impossible to differentiate one from 
the other by eye. What is different is shown as thin lines corresponding to the individual contributions of the Hi and Dl ISM components. 
Vertical tick marks give the spectral positions of each component. The reference stellar continuum Sr, (thick dotted line) as evaluated by 
Wood (2001) from the adjustment of several data sets is shown along with the stellar continuum (Sr corrected by a 8 th order polynomial, 
thick solid line) found in the fitting process. If the stellar profile is not known to better than ±10%, the Hi column density on the line of 
sight is poorly evaluated as well as the D/H ratio. 



